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Abstract 

The viability of neutralino dark matter as an explanation of the 130 GeV gamma ray 
signal from the Galactic Center recently observed by the Fermi Large Area Telescope is 
examined. It is found that the signal can be compatible with a sharp feature from inter- 
nal bremsstrahlung from a mostly bino dark matter particle of mass around 145 GeV, 
augmented by a contribution from annihilation into gamma+Z via a small wino ad- 
mixture. This scenario circumvents the problematic overproduction of lower energy 
continuum photons that plague line interpretations of this signal. Sleptons approxi- 
mately degenerate in mass with the neutralino are required to enhance the internal 
bremsstrahlung feature. 



1 Motivation 



Several analyses ( PQ [2] [S] |1] ) have recently confirmed the presence of a sharp feature, in- 
compatible with conventional astrophysics, at an energy of approximately 130 GeV in the 
gamma ray spectrum towards the Galactic Center in the data gathered by the Fermi Large 
Area Telescope (LAT). While the possibility that this might be an instrumental effect or a 
product of nonconventional astrophysics still exists, such a feature has long been earmarked 
as a "smoking gun" signature of dark matter annihilation in the galaxy; this tantalizing 
interpretation has therefore generated significant excitement. 

Assuming a dark matter origin, the signal is best fit by a 130 GeV dark matter par- 
ticle pair-annihilating into photons with an annihilation cross section of (crv)^^ = 1.27 x 
10~^^cm^s~^, assuming an Einasto profile for the dark matter distribution |T]. From a par- 
ticle physics point of view, this scenario poses two major problems. First, since dark matter 
is not expected to couple directly to photons, annihilation to a photon pair must occur via a 
loop (see [5],[6j for a full calculation of this process in supersymmetry) ; for a thermal relic, 
this loop-suppressed cross section is generally too small to produce the signal observed by 
Fermi. Second, even if this cross section can be made large enough, tree-level annihilation 
to particles that mediate the photon pair production process should produce a large contin- 
uum of photons at lower energies, which is not seen in the Fermi data [7J[8J[9J[10J. These 
considerations have been shown to rule out the most promising and the most studied dark 
matter candidate, the lightest neutralino in supersymmetry, as an explanation of this line 
signal 

However, a monochromatic line signal is not the only possibility that can explain this 
feature; it is well-known that internal bremsstrahlung (hereafter IB) - the production of a 
photon in conjunction with the leading annihilation channel into charged particles - can also 
give sharp spectral features in the 7 ray spectrum close to the dark matter mass [TT] [12] [13] . 
The first evidence of the 130 GeV feature was, in fact, found in searches for IB signals in 
the Fermi LAT data [1], where a fit significance of 3.2(4.6)(T excess with (without) the look 
elsewhere effect was found for an IB signal corresponding to annihilation of a dark matter 
particle of mass around 150 GeV. While an IB signal is broader and gives a poorer fit to 
Fermi data than a line signal, it already appears at tree-level and can potentially overcome 
the constraints from signal strength and overproduction of continuum photons that plague 
line signals. 

The purpose of this paper is to examine whether IB from neutralino dark matter an- 
nihilation can provide a viable explanation of the 130 GeV signal observed by Fermi. The 
possibility of explaining the feature with IB has been mentioned in several papers, but an 
explicit treatment within supersymmetry is still lacking. While the toy model employed 
in [4j is very similar in several respects, there are several crucial differences, such as the 
availability of several IB channels and significant contributions from the 77 and 7Z lines, 
and the supersymmetric parameter space is not treated in detail. Given the overwhelming 
popularity of the neutralino as a dark matter candidate, a more careful study is therefore 
desirable. 
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This paper is organized as follows. Section 2 discusses continuum constraints and IB in 
the context of supersymmetry. Section 3 contains details and discussions of the scans, a few 
benchmark points, and their fit to the Fermi data. Section 4 discusses relic density and other 
relevant dark matter considerations. The main results are summarized in Section 5. 

2 Neutralino Dark Matter and Internal Bremsstrahlung 

In the Minimal Supersymmetric Standard Model (MSSM), neutralino dark matter is an 
admixture of the bino, neutral wino, and the two neutral Higgsinos, and its composition 
determines the strength of the gamma ray signal and the relative size of the continuum. 
Table [T] lists the approximate annihilation cross sections into 77, 7Z, and the total cross 
section for a dark matter particle of mass 130 GeV in the three limits (taken from |7j). 



Model 


a^^v (cm^s ^) 


(J^zv (cm^s ^) 


(TtotalV (cm^S 1) 




Bino 

Wino 

Higgsino 


~ 10-30 

2.5x10-2^ 
1.1x10-28 


~ 10-31 
1.4x10-26 
3.7x10-28 


~ 10-27 

4x10-24 
4.2x10-25 


~ 1000 
210 
710 



Table 1: Annihilation cross sections for various neutralino constituents. The bino cross 
sections are with slepton masses at 200 GeV. R^^ represents the size of the continuum relative 
to the line signal [7|. 

Recall that the best fit to the Fermi data requires an annihilation cross section of a^^v = 
1.27 X 10-27 cm^s-i for an Einasto profile. While the wino and Higgsino line signals are at 
the right order of magnitude to produce this signal, the bino line signals fall several orders of 
magnitude short of this requirement. The required cross section can be lowered by allowing 
for a steeper dark matter profile at the Galactic Center - for instance, |1] finds that a cross 
section of ~ 2 x 10-2^ cm^s-^ can explain the signal close to the Galactic Center (Reg4 and 
Reg5 in [1]), where the signal is the most significant, for a contracted Navarro, Frenk, and 
White (NFW) dark matter profile with slope a = 1.3. This, however, is still insufficient to 
bring line signals from the bino into contention. 

The wino and Higgsino, meanwhile, annihilate dominantly into gauge bosons, and their 
subsequent decays produce a significant continuum of photons at lower energies. The size of 
this continuum is represented by the ratio i?*'^, listed in the final column of Table [l| defined 
in [7] as 

j^th _ ^onn ^ 

which is to be constrained through comparison with the analogous ratio from observation [7j 



where N refers to the number of photons from the relevant process, and is the total 
number of photons per annihilation in the considered energy range. In [7], simply requiring 
that the continuum contribution not supersaturate the data was found to constrain to be 
below 75 to 12C[^for dark matter in the mass range 125 to 150 GeV annihilating primarily into 
W or Z bosons. Line contributions from wino or Higgsino dark matter are in clear tension 
with this bound, ruling them out as an explanation of the feature observed by Fermi. It 
should be noted that a 130 GeV dark matter particle annihilating into Ws is also in tension 
with PAMELA antiproton constraints (see e.g. and observation of dwarf galaxies 

([IH]; also see y/Tj for related uncertainties with dwarf galaxies). 

Despite the incompatibility of a line signal, internal bremsstrahlung (IB) might, as moti- 
vated in the previous section, still salvage supersymmetry as an explanation of the 130 GeV 
signal observed by Fermi. IB refers to radiation of a photon from either the charged Standard 
Model (SM) final states that dark matter annihilates into or the charged mediator in the t- or 
u-channel. One can nominally distinguish between the two as final state radiation (FSR) and 
virtual internal bremsstrahlung (VIB) respectively, but they cannot be treated separately in 
a gauge-invariant manner and must always be considered and calculated together. 

The IB component from neutralino annihilation is known to be the most prominent when 
annihilation is into particles that are effectively massless relative to the neutralino, and the 
virtual particle that mediates the process is close in mass to the neutralino [II][I2]- Since 
the W and Z gauge bosons are massive final states for a neutralino of mass around 130- 
150 GeV, IB from wino or Higgsino dark matter does not produce a feature sharp enough to 
explain the Fermi observation, despite the presence of a degenerate chargino to mediate its 
annihilation; this has been verified explicitly. 

IB from bino dark matter, on the other hand, is more promising. The main annihilation 
channels for a bino are to fermion pairs, mediated by the corresponding sfermions. For 
nonrelativistic annihilation in the halo, the cross section for this process is helicity suppressed 
by a factor of {rrif /m^)'^] since the top quark is heavier than the dark matter mass of interest 
here and all other SM fermions are 0{GeV) or lighter, this suppression is of several orders of 
magnitude, and acts as an efficient mechanism to suppress the continuum photon production. 
The addition of a photon in the final state, on the other hand, lifts this helicity suppression, 
and crv{xx ~^ f fl) can be comparable to avixx ~^ ff)- Since fermions and sfermions 
couple to the bino via hypercharge and leptons have larger hypercharge than quarks (also, 
sleptons are generally lighter than squarks), IB primarily involves leptonic channels. 

For an almost pure bino, the IB cross section is fairly robust; in the limit of massless 
fermions, it is approximately given by [llj 



da. 



dx 



X 




(1 



{l + ^i-2x) 



x) 



(1 +/i)(l +/i-2x) 
{l + ^i-xf 



log 



1 + li-2x 



1 + ^l 




Allowing for a power law background makes this constraint even stronger [7]. 
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where /x = {m^-^/m^y = {mj-^/m^'^ (assuming the same mass for both sfermions), gL{R) 
corresponds to the couphng of the left(right) handed sfermion to the bino, and x = E^/m^. 
Note the prominent log enhancement close to the kinematic edge {x ~ 1) when sfermions 
are approximately degenerate with the bino (/i ~ 1); light sleptons are therefore a crucial 
element of strong IB features. 

This setup offers a clear strategy towards an attempt at a supersymmetric explanation 
of the Fermi signal. The primary contribution must come from IB from mostly bino dark 
matter in the 130 to 150 GeV range, with sleptons not too far above in mass to produce a 
sufficiently large and sharp feature. A small wino and/or higgsino component can augment 
this signal via a 77 or 7Z peak as long as the production of continuum photons is sufficiently 
suppressed to evade the observational bounds. 

3 Approach, Results, and Discussion 

The following analysis is based on Fermi LAT data from the inner 3° radius region around 
the Galactic Center, where the 130 GeV signal was found to be the most significant |2J [3J . For 
this purpose, event counts as listed in Appendix A of (unmasked region), corresponding 
to ULTRACLEAN events in the Pass 7_Version 6 releasfl for 128 energy bins from 5.1 GeV 
to 198 GeV are used. The interested reader is referred to that paper for details of energy 
binning and individual photons counts. 

The gamma ray spectrum from neutralino annihilation is generated using DarkSUSY 
version 5.0.5|^[T8]. The spectrum is normalized such that a^^v = 2x 10~^^cm^s~^ corresponds 
to 30.3 photons in the 77 peak. This choice is made for the following reason. The authors of 
[7] find that the data set in question is best fit by a 77 peak with 30.3 photons. The lowest 
cross section into two photons that is consistent with the signal in this region, meanwhile, 
is roughly a^^v = 2 x lO^^^cm^s"^ (^Lj, for Reg4 and Reg5 therein), for a contracted NFW 
dark matter profile with a modified slope a = 1.3. While an Einasto profile with a^^f = 
1.27 X 10~^'^cm^s~^ is more consistent with data from extended regions beyond the 3° region 
considered here [Ij the number chosen above is a better refiection of what the cross section 
at least needs to be in order to be considered a viable option. 

Corrections for the Instrument Response Function of the Fermi LAT, which describes the 
energy dispersion of incident photons, are crucial for line signals and the sharp IB feature 

^http: / /heasarc. gsfc.nasa.gov/FTP / fcrmi /data/lat / weekly / p7v6 
A subtlety regarding the IB component from DarkSUSY is worth mentioning here. DarkSUSY obtains 
the spectrum from Pythia, which simulates the annihilation process as a decay of a hypothetical particle of 
mass 2m^, hence missing contributions corresponding to virtual IB. DarkSUSY corrects for this by subtract- 
ing the FSR spectrum of the hypothetical decay from the Pythia result and adding the full IB contribution 
[TT] . Spectra obtained from Pythia and DarkSUSY might show differences due to this correction. The author 
thanks Joakim Edsjo for clarification on this issue. 

""^See [3D] for a detailed discussion of the compatibility of various profiles to the signal from extended 
regions. In particular, the NFW profile used here is consistent with the signal in the 3° region considered 
here but needs to be modified at larger angles to maintain consistency with data from extended regions. The 
author thanks Torsten Bringmann and Christoph Weniger for pointing this out. 
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at the kinematic edge, and less important for the continuum at lower energies. This is 
taken into account by fitting the energy dispersion at 130 GeV, plotted in Appendix B of |[7] , 
with a Gaussian, and applying this dispersion correction to the dark matter spectrum above 
100 GeV. Correction for the change in effective area of the instrument at different energies 
is also approximately incorporated using the information provided in [19]. A more careful 
treatment of these factors, while possible, is unnecessary for the major objectives of this 
paper. 

Finally, in addition to the continuum from the decay of annihilation products, other 
potentially important dark matter contributions to the gamma ray spectrum also need to 
be considered. Of these. Inverse Compton Scattering (ICS) of photons in the interstellar 
radiation field off charged products from dark matter annihilation is the most important. 
The ICS contribution is estimated using a semi-analytic formalism described in [20].|21j. 
with the simplification of ignoring spatial diffusion (equivalent to setting the halo function 
to unity), and found to be negligible; this is understandable, since the production of fermions 
is strongly helicity suppressed. 

3.1 Scan Results 

A scan over the MSSM parameter space, optimized for bino dark matter and the Fermi 
signal, was performed with DarkSUSY. The lightest neutralino was required to have a mass 
between 120 GeV and 160 GeV, and be mostly bino, with the wino/Higgsino components 
required to be sufficiently suppressed to avoid overproduction of the continuum. A single 
mass value was chosen for all sleptons, and this value was constrained to be within 20 GeV 
of the lightest neutralino mass. Squark masses were fixed at 1 TeV. The bino mass Mi was 
varied between 120 and 160 GeV, and the wino mass M2 and the higgsino mass parameter 
fi were varied between 120 GeV and 2 TeV; all three were allowed to take negative values. 

In the Fermi data set used in this analysis, the excess appears at energies between 
121.62 GeV and 136.40 GeV, where a total of 24 photons are observed (Appendix A in [7]). 
The first issue of concern is whether the signal from neutralino annihilation is strong enough 
to explain this excess. A falling power law background, obtained from a fit to the spectrum 
at lower energies, contributes 6 or 7 photons in this energy range. A good fit should be pos- 
sible with ~9 or more photons from neutralino annihilation in addition to this background 
contribution, as the observed count of 24 is then 2a or less away (assuming cr ~ V^)- For 
comparison, a 77 peak with 30.3 photons, the best monochromatic fit to data [7J, contributes 
~19 photons in this energy range after energy dispersion. 

The first plot in Figure [T] shows the number of photons in this energy range between 
121.62 GeV and 136.40 GeV from dark matter annihilation as a function of dark matter mass. 
The number is fairly robust in the range of dark matter masses that can explain the Fermi 
signal: dark matter contributes 0{few) photons, peaking at ~ 13 photons around ~ 
145 GeV. The occurrence of the peak at this energy is understandable: for ~ 145 GeV, 
the 7Z channel gives monoenergetic photons at Ej = m^{l — m\ 

/Ami) ~ 130 GeV, which 

can be sizable even with a tiny wino/Higgsino component. Meanwhile, the peak of the IB 
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Figure 1: The total number of photons in the 121.62 GeV to 136.40 GeV energy range from 
dark matter annihilation, as a function of dark matter mass (left) and slepton-neutralino 
mass difference (right), for a scan over primarily bino candidates. Points with thermal relic 
density calculated to be in the range 0.1 < Qh"^ < 0.124, making them consistent with relic 
density constraints, are shown in black. 



feature still falls mostly into the 121.62 GeV to 136.40 GeV range (recall that the best fit for 
a purely IB signal occurs for ~ 150 GeV); hence both IB and the jZ line are at ideal 
energies to contribute towards the signal. As the mass changes away from this ideal value, 
either the IB peak or the line contribution is lost, and the signal dies away, as is seen in the 
plot. 

The second plot in Figure [T] shows the dependence on the slepton-neutralino mass differ- 
ence, which follows what is expected from Equation [3] the photon count gradually rises as 
the sleptons become more degenerate in mass with the neutralino. The scan is terminated 
at 20 GeV, as the downward trend continues and a sufficient number of photons cannot be 
obtained. 

The results of this scan show that, with a combination of sleptons approximately degen- 
erate with the dark matter particle and contributions from both IB and line components, a 
dark matter signal large enough to explain the Fermi feature can be possible without over- 
production of the continuum. Attention must now switch to whether the contribution is of 
the right shape to explain the observed signal. 

3.2 Benchmark Points and Fit to Data 

This section discusses four benchmark points (labelled BMl, BM2, BM3, and BM4) that 
are representative of the scanned sample, and their fits to the Fermi signal. These points 
are listed in Table |2| with other relevant information. It should be stressed that these were 
chosen to highlight distinct features of signals that are possible with internal bremsstrahlung, 
and are not the points that best fit the data. 

To perform the fit to Fermi data, the spectrum generated with DarkSUSY for each 
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BMl 


BM2 


BM3 


BMi 


Ml 


135.2 


144.7 


145.6 


138.2 


M2 


235.5 


152.8 


150.4 


161.2 




-489.9 


838.4 


783.0 


512.9 


tan/3 


18.5 


6.6 


33.2 


20.5 




136.7 


156.6 


146.7 


138.5 




134.4 


143.0 


144.7 


136.4 


bino fraction 


0.99 


.90 


0.91 


0.97 




0.19 


.0058 


0.0033 


0.11 


ra^ from IB 


4.8 


1.8 


4.5 


14.7 


(77 + iZ) 


2.0 


5.1 


5.2 


4.4 


TS 


15.8 






17.8 


significance 


4.0 






4.2 



Table 2: Four benclimark points cfiosen for detailed study and fit to Fermi data, and fit 
results. All masses are in GeV. refers to the number of photons contributed to the 
121.62 GeV to 136.40 GeV energy bin. 



benchmark point is added to a single falling power law background, with the normalization 
and the power law index allowed to vary to give the best fit. Such fits are not possible for 
points for which the low energy continuum is close to saturating the Fermi data; however, 
such points are still useful to illustrate various characteristics of IB signals. This is true 
of benchmark points BM2 and BM3. Therefore, in this section, fits with background are 
only performed where feasible (BMl and BM4). It should be noted, however, that the 
assumption of a single power law background across the entire spectrum is a rather strong 
one, and signals not allowed by such a background might be consistent with other forms of 
background. 

Following |l],[Tj, the significance of the fit is estimated by maximizing the likelihood 
function 

^bins 

In £ = ^ nfc ■ In (pk - (pk - ln(nfc!) (4) 

k=l 

where (0^) represents the number of photons observed (expected) in the k^^ energy bin. 
This is then used to calculate the test statistics (TS) 

TS =-2ln (5) 

•'-benchmark 

where /^benchmark IS the likelihood of the benchmark model, and £nuii corresponds to the null 
hypothesis, i.e. fit with a power law background only. The nominal significance of the fit is 
taken to be y/TS. 

The gamma ray spectra from each of the benchmark points, superimposed on the Fermi 
data over the entire energy range between 5 and 198 GeV, are plotted in Figure [2j Informa- 
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Figure 2: Fermi data from the inner 3° of the Galactic Center for all 128 energy bins (blue 
dots, as listed in Appendix A of [7]) and the gamma ray spectra from dark matter (green) 
for the four benchmark points BMl, BM2 (top row), BM3, and BM4 (bottom row). The 
black curves for BMl and BM4 represent the overall fit to the signal, consisting of a single 
power law background in addition to the dark matter signal; BM2 and BM3 supersaturate 
the continuum at lower energies and do not allow such fits. 

tion regarding the significance of the fit is listed in Table |2j As mentioned earlier, the fit is 
only performed for two of the benchmark points, BMl and BM4; the other two benchmark 
points contain large continuum signals at low energies that saturate the Fermi data and are 
therefore incompatible with a single power law background, making such fits impossible. 
The reader is advised to use caution in interpreting these fit results, since these are not pro- 
duced from an extensive scan and detailed fitting procedure, which is not the main purpose 
of this paper, and are merely meant to be a rough indication of the compatibility between 
prediction and signal. 

Next, each benchmark point and its fit to data is discussed in turn. 

BMl is an almost pure bino that contributes dominantly via IB. The mass degeneracy 
between the neutralino and the sleptons - the difference is only 2 GeV - makes the IB 
feature very sharp and prominent above the rest of the spectrum. This very mass degeneracy 
also facilitates coannihlations in the early universe, leading to a relic density very close to 
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the observed value. It can be seen that the pure bino comes with an extremely efficient 
suppression of continuum photons - the continuum is so small that it does not even appear 
on Figure |2] - due to helicity suppression of annihilation into fermions. However, it also 
comes with the disadvantage that the desired number of photons is difficult to obtain via a 
purely IB contribution, even with extremely degenerate sleptons. 

BM2 improves on BMl by introducing a 10% wino admixture to the neutralino, opening 
up significant contributions from the 7Z line: in this case, the line contributes 5 photons 
to the signal bin. On the other hand, this wino contribution also floods the spectrum with 
continuum photons at lower energies that, for this particular benchmark point, saturate the 
Fermi signal in the region around 15—20 GeV (Figure[2]), making the signal incompatible with 
a single power law background. The thermal relic density also plummets due to extremely 
efficient annihilation through the wino component. The sleptons are about 15 GeV heavier 
than the neutralino, and the lifting of this mass degeneracy suppresses the IB contribution 
relative to BMl, so that the line signal is the primary component contributing to the 130 GeV 
signal. Although BMl and BM2 contain essentially the same number of photons in the signal 
bin, the line provides a narrower feature that better fits the shape of the Fermi signal around 
130 GeV. 

BM3 combines the virtues of both BMl and BM2: sleptons degenerate with the neu- 
tralino lead to a sharp IB feature, while a small wino admixture contributes a prominent 
7Z line signal, resulting in a relative abundance of photons in the right energy bin. Such 
points, characterized by sleptons approximately degenerate with the neutralino as well as 
line contribution from a small wino component, should be considered the most promising 
avenues for producing the observed signal. As with BM2, dark matter contribution saturates 
the Fermi signal in the 15 — 20 GeV window, making the signal incompatible with a single 
power law background, and the thermal relic density is smaller than desired due to the wino 
component. 

BM4 analyses a fit with about 19 photons - the same number as from a best fit monochro- 
matic line - in the 121.62 GeV to 136.40 GeV energy bin. This is achieved by enhancing the 
dark matter signal by an additional factor of 3, and could correspond to some astrophysical 
enhancement such as a steeper dark matter profile at the Galactic Center or substructure 
along the line of sight towards this region. A significantly better fit is obtained in this 
case, and the enhancement allows the signal to be composed mostly of IB, enabling further 
suppression of the wino component relative to BM2 or BM3. This suppression has two ad- 
vantages: one, the continuum is suppressed relative to BM2 and BM3, allowing for a good fit 
with a single power law background; two, the relic density is raised to a value that matches 
observation. 

That the subdominant admixtures in all of BM2, BM3, and BM4 are winos is not a 
coincidence. There are several reasons for the wino being preferred to Higgsinos. Because 
Higgsinos fare far worse with continuum photons than winos (see Table |l|), they can con- 
tribute no more than 1 or 2 photons to the 130 GeV signal before saturating the spectrum at 
lower energies. In addition, since the Higgsino annihilation cross sections - both total and 
into lines - are smaller than the corresponding ones for a wino by an order of magnitude. 
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producing the same number of line photons as the wino requires a larger Higgsino fraction, 
suppressing the bino fraction and consequently the IB contribution. Finally, the 7Z line is 
more prominent than 77 for a wino, which is desirable. 

The four benchmark points discussed in this subsection each illustrate an important 
advantage - suppression of continuum photons, elimination of the need for sleptons within 
a few GeV of the lightest neutralino, maximization of photon count in the bins where Fermi 
observes a signal, and significant improvements from additional 0{1) boosts, respectively - 
but demand caution in other aspects. The above analysis serves to highlight the interplay of 
the various factors that are relevant for a supersymmetric explanation of the Fermi signal, 
but confirms nevertheless that such an explanation is indeed possible with IB. 

4 Other Considerations 

This section is devoted to brief discussions of various related aspects that deserve attention. 

4.1 Relic Density 

The bino, being a gauge singlet, is generally inert, leading to a thermal relic density that, 
for a bino of mass around 130 GeV, is far in excess of the observed value. However, there 
exist well-understood ways to resolve this discrepancy. One is to have an almost degenerate 
slepton (usually stau), within ~ 5% of the dark matter mass, to enable efficient coannihilation 
|22j : another is to make the neutralino well-tempered, i.e., introduce a small admixture of 
wino or Higgsino [2S]- Both of these features are prominent in the class of models discussed 
in this paper (recall that the former led to a reasonably good relic density for BMl, while 
both factors were at play in obtaining the correct relic density for BM4), and obtaining a 
thermal relic density in agreement with observations appears to be a tractable task. For 
the set of points scanned over with DarkSUSY, the computed relic densities ranged from 
0.001 < Qh"^ < 0.3, with a significant fraction within 2a of the current best fit value; recall 
that these points that are in agreement with observation were plotted in black in Figure [T] 

4.2 Collider and Direct Detection Constraints 

Collider and direct detection searches generally place stringent constraints on dark matter 
and supersymmetry, and must be considered. For mostly bino dark matter and squarks 
at the TeV scale or heavier, Tevatron and LHC constraints are easily avoided. Despite 
the candidate considered here having appreciable couplings to leptons, facilitated by the 
presence of light sleptons, LEP constraints on dark matter [24j only apply to dark matter 
masses below its threshold of 100 GeV, and are irrelevant to the mass range of interest in this 
paper. Likewise, when squarks are heavy and the lightest neutralino is a sufficiently pure 
gaugino, tree level spin-independent direct detection interactions with nuclei are suppressed, 
and the candidate is safe from the direct detection bounds placed by XENONIOO [25j (see 
[26], [27] for elaboration and a detailed study of this point). The dark matter candidates 
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studied in this paper are therefore safe from both coUider and direct detection constraints 
at present. 

4.3 Astrophysical Uncertainties 

While gamma rays represent the cleanest indirect detection channels for dark matter, a 
dark matter interpretation of the Fermi 130 GeV signal is still plagued with astrophysical 
uncertainties. Of these, a precise knowledge of the astrophysical background - even at the 
level of whether it follows a single power law across the entire energy range of interest - 
and reliable knowledge of the dark matter profile at the Galactic Center, or the existence 
of substructures in the direction of the signal, introduce huge uncertainties in translating 
observations to implications for possible underlying particle physics models. For instance, 
if dark matter at the Galactic Center follows an isothermal profile, the combination of IB 
and line signals presented here are no longer plausible as an explanation of the Fermi signal. 
On the other hand, even a modest presence of dark matter substructure in the direction of 
the inner 3° of the Galactic Center can greatly enhance a dark matter signal, allowing a 
larger portion of supersymmetric parameter space - such as purely IB contributions, or IB 
with heavier sleptons - into consideration, or enables better fit^ even with less peaked dark 
matter profiles such as the more favored Einasto profile. It should therefore always be kept 
in mind that the uncertainties introduced by astrophysical factors are large and can have 
significant implications. 

4.4 Naturalness 

With a relatively heavy Higgs discovered and the LHC failing to find any light superpartners, 
the naturalness of supersymmetry as a resolution of the hierarchy problem has become an 
important issue. A commonly employed measure of fine-tuning in terms of the tree level 
Z boson mass requires fi ~ rriz, with the amount of fine-tuning scaling as ~ {fi/rrizY (see 
e.g.[2H],[2S] for more detailed discussions). The scenario presented in this paper, where the 
neutralino is almost entirely gaugino, and the /j, parameter is required to be extremely large 
in order to suppress continuum contributions from the Higgsino component, is significantly 
fine-tuned in this regard. This, however, is strictly true only in the MSSM, where the 
fine-tuning problem is already known to be serious. In contrast, in nonminimal versions of 
supersymmetry favored by a 125 GeV Higgs[29j, such as the Next-to-Minimal Supersymmet- 
ric Standard Model (NMSSM) or A-SUSY, a parametric suppression of fine-tuning can occur, 
and larger values of /i can be perfectly natural (see [27] for a detailed discussion). Although 
the discussion presented here was confined to the MSSM, a neutralino that is mostly bino 
with a small wino component can be easily realized in such nonminimal extensions, and 
would be consistent with naturalness. 

^BM4 serves as an illustration of this point. 



11 



5 Conclusions 



There is now clear evidence of an unexplained feature at 130 GeV towards the Galactic Center 
in the Fermi LAT data. A dark matter interpretation is very tempting, and the feature fits 
extremely well to dark matter annihilating into monochromatic photons, the long anticipated 
"smoking gun" signature of dark matter. The most studied dark matter candidate, the 
lightest neutralino in supersymmetry, is incompatible with this line interpretation of the 
signal, constrained by the absence of a continuum at low energies in the observed data, or 
small cross sections into line photons. 

The purpose of this paper was to examine whether these constraints can be circum- 
vented in supersymmetric scenarios where internal bremsstrahlung plays a prominent role. 
Satisfactory scenarios were indeed found. A few benchmark models illustrating the major 
possibilities were presented, and their agreement with Fermi data explored. 

The scenario most consistent with the 130 GeV signal corresponds to sharp internal 
bremsstrahlung from a ~ 145 GeV mostly bino dark matter particle in conjunction with 
a 7Z line from a subdominant wino component; with this choice of mass, the peaks of both 
the IB and 7Z signals fall in the 130 GeV region, producing a strong signal. Light sleptons 
approximately degenerate with the neutralino are required to make the IB feature promi- 
nent. This combination of bino dominance, approximate mass degeneracy of the sleptons 
and the neutralino, and a possible line contribution from a subdominant wino component is 
a generic feature of the class of candidates studied in this paper. The presence of light slep- 
tons also facilitates coannihilations, providing thermal relic densities roughly in agreement 
with observation. 

A contracted NFW profile was chosen over the more favored Einasto profile to allow for a 
more generous - therefore broader - treatment of the parameter space; the photon flux from 
choosing the latter profile is only an (9(1) factor smaller, and a modest contribution from, for 
instance, substructure along the line of sight to the center of the galaxy can easily overcome 
this difference. Given the large uncertainties in these astrophysical factors, possibilities of 
such corrections should not be ignored. 

While instrumental or nonconventional astrophysical effects might yet explain this 130 GeV 
anomaly, the possibility that this might be the first signature of dark matter - of a particle 
beyond the Standard Model - is one with tremendous implications, and one worth pursuing 
even in the midst of uncertainty. More data, from the Galactic Center and elsewhere, and 
with Fermi as well as with other instruments, will gradually improve the details of the signal, 
leading to a clearer picture. For the moment, the possibility that the signal has its origins 
in dark matter annihilation remains alive; this paper has presented a case that so too does 
the possibility that that origin is supersymmetric. 

Note Added 

During the completion of this project, a review paper [SD] appeared, where an MSSM scan 
exploring the relative sizes of the IB, line signals, and secondary photons in the context of 
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the Fermi 130 GeV signal is presented and discussed. The resuhs presented in this paper are 
in agreement with the results therein. 
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